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ABSTRACT 

We use a one-dimensional hydrodynamical code to investigate the effects of preheating 
on gas accretion and cooling in cold dark matter halos. In the absence of radiative 
cooling, preheating reduces the amount of gas that can be accreted into a halo, and 
the accreted gas fraction is determined by the ratio of the initial specific entropy of 
the gas to the virial entropy of the halo, S p h/S v . In the presence of radiative cooling, 
preheating affects the gas fraction that can cool in two different ways. For small halos 
with masses M < 10 12 h~ x M Q , preheating suppresses gas accretion, but most of the 
accreted gas can cool. For more massive halos, preheating affects the cold gas fraction 
not only by reducing the amount of accreted gas, but also by reducing the cooling 
efficiency. For both small and massive halos, gas cooling is delayed by preheating 
if the halo gas is assumed to be a single-phase medium. However, if the halo gas is 
assumed to be a multi-phase medium, cooling can occur over a wider range of redshifts. 
Unlike in a single-phase medium where cooling is inside-out, gas cooling in a preheated 
multi-phase medium can occur simultaneously over a wide range of radii. As examples, 
two specific preheating cases are investigated. In the first case, the preheating specific 
entropy is assumed to be proportional to the virial entropy of the halo, Sph oc S v , as 
expected from AGN feedback. Such preheating effectively suppresses radiative cooling 
in halos with M > 10 13 h^ 1 M , but has little effect on smaller halos. We suggest 
that this may be the reason why the stellar mass function of galaxies breaks sharply 
at the massive end. Such preheating also helps create the hot diffused halos within 
which the "radio mode" feedback of AGNs can act effectively. In the second case, the 
intcrgalactic medium is assumed to be warm. Here the total amount of gas that can 
cool in a halo scales with halo mass as oc M 2 , as would be required to match the 
observed stellar- and Hi-mass functions in the current CDM model at the small mass 
end. 
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1 INTRODUCTION 

According to the current theory of structure formation, 
galaxies are believed to form in cold dark matter (CDM) 
halos that are virialized clumps of CDM formed through 
gravitational instability from initial density perturbations 
in the early universe. In this scenario, the problem of galaxy 
formation can be divided roughly into two parts: (1) the 
formation of CDM halos in the cosmic density field, and 
(2) the accretion and cooling of gas into dark matter halos 
to form stars and the structures of luminous galaxies. The 
buildup of dark matter halos is now fairly well understood 
using both analytic models (e.g. Press & Schechter 1974; 
Bond et al. 1991; Sheth, Mo & Tormen 2001) and A-body 
simulations (e.g. Navarro, Frenk & White 1996; Jing & Suto 
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2000; Moore et al. 1998; Zhao et al. 2003a, b; Wechsler et 
al. 2002) . However, how gas cools and assembles to form the 
luminous part of a galaxy is still an open question. 

The pioneering work by Binney (1977), Rees & Ostriker 
(1977) and Silk (1977) assumed that the gas in a dark mat- 
ter halo is heated by shocks as it collapses into dark matter 
halos. The heated gas then cools via radiative processes, and 
forms stars in the centers of dark halos. This assumption has 
been combined with the halo population predicted by CDM 
models to make predictions for the properties of the galaxy 
population, such as the luminosity function, the colour dis- 
tribution, etc (White & Rees 1978; White & Frenk 1991). 
With modern computational facilities and numerical codes, 
many of the processes in galaxy formation can be studied in 
more detail with numerical simulations (e.g., Katz & Gunn 
1991; Cen & Ostriker 1992; Springel & Hernquist 2003). 

Despite the progress made so far, galaxy formation re- 
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mains one of the most challenging problems in astrophysics. 
One of the long-standing problems in the CDM scenario of 
galaxy formation is to match the halo mass function with the 
luminosity function of galaxies. The galaxy luminosity func- 
tion has a characteristic shape that is very different from 
the halo mass function. The faint end slope of the galaxy 
luminosity function is significantly shallower than the halo 
mass function at the low-mass end. Also, the galaxy lumi- 
nosity function has an abrupt break at the bright end, while 
such decline in the halo mass function occurs at a mass that 
corresponds to a much lower abundance (Yang et al. 2003). 
Furthermore, as shown in Yang et al. (2005), the relation 
between the luminosity of the central galaxy in a halo and 
the halo mass is strongly nonlinear at both the low-mass 
and high-mass ends. These results imply that the efficiency 
of star formation must be suppressed in both low-mass and 
high-mass halos. 

Various suggestions have been made to explain the de- 
pendence of star formation efficiency on halo mass. For mas- 
sive halos, it has been suggested that the low star formation 
efficiency is a result of the inefficient gas cooling. Thoul & 
Weinberg (1995) studied gas cooling in massive halos in con- 
siderable detail and found that a significant amount of gas 
can still cool in such halos, although the cooling efficiency is 
lower than that in lower mass halos. Thus, radiative cooling 
alone may not be able to explain the bright-end behavior in 
the galaxy luminosity function. For low mass halos, sugges- 
tions have been made that heating by the UV background 
and/or supernova explosions associated with star formation 
may suppress gas cooling and hence the star formation effi- 
ciency. However, heating by the UV background is effective 
only in small halos with masses M < 10 h~ l Mq at low 
redshift (Quinn et al. 1996; Gnedin et al. 2000; Hoeft et 
al. 2006) and its effect is reduced by self-shielding at high 
redshifts (Dijkstra et al. 2004). Heating by supernova explo- 
sions may also be ineffective because of the development of 
Ray leigh- Taylor instabilities in the interstellar medium (e.g. 
Mac Low & Fcrrara 1999; Strickland & Stevens 2000). 

One important process that can prevent gas from cool- 
ing too fast in dark matter halos is preheating (e.g. Mo & 
Mao 2002; Oh & Benson 2003; Scannapieco & Oh 2004; Mo 
et al. 2005). In this scenario, the intergalactic gas is assumed 
to be heated to some finite entropy before it is accreted into 
dark matter halos. This will affect how gas is accreted into 
dark matter halos (e.g. Mo & Mao 2002; van den Bosch, 
Abel & Hernquist 2003) , and the distribution, thermal state 
and radiative cooling of the accreted gas. In this paper, we 
use 1-D simulations to study in detail the effects of pre- 
heating on galaxy formation in dark matter halos. We find 
a number of scaling relations that are useful to understand 
and to model galaxy formation in preheated media. Further- 
more, we consider two realistic cases of preheating and study 
their consequences on the formation of the central galaxies 
in dark matter halos. In one, preheating is assumed to be 
generated by the energy feedback of AGNs. We show that 
such preheating can effectively suppress gas cooling in halos 
more massive than ~ 10 13 M . In the second case, we con- 
sider a constant, relatively low entropy floor and study its 
effect on mass accretion and cooling in low-mass halos. This 
model is motivated by the fact that low-mass halos at low 
redshift may be surrounded by a warm intergalactic medium 
generated by gravitational pancaking (Mo et al. 2005). Our 



results show that the cold gas fraction in small dark matter 
halos can be significantly reduced in such a warm medium. 

This paper is organized as follows. We briefly describe 
the numerical scheme of our 1-D hydrodynamical code and 
simulation models in Section 2. In Section 3, we show adia- 
batic simulation results and discuss the effects of preheating 
on gas accretion. In Section 4, we study the cooling of pre- 
heated media in halos of different masses using both single- 
phase and multi-phase simulations. We implement our sim- 
ulations to study the two particular preheating models in 
Section 5. Further discussion and a summary of our results 
are given in Section 6. 



2 MODEL AND SIMULATION 

In this paper, we use 1-Dimensional Lagrangian hydrody- 
namical simulations to study the accretion and cooling of 
gas in dark matter halos. In these simulations, the mass dis- 
tribution is assumed to be spherically symmetric, and the ra- 
dial distributions of the dark matter and gas are represented 
by spherical shells, each of which contains a certain amount 
of mass. The motions of the dark matter shells are followed 
with the same 1-d code as described in Lu et al. (2006). In 
what follows, we describe how we simulate the evolution of 
the gas component. 



2.1 Hydrodynamics and thermodynamics 

We treat the gas component as fluid, so the motion is gov- 
erned by the following set of equations: 
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These are the continuity equation, the momentum equation, 
the energy equation, and the equation of state, respectively. 
In these equations, p g , p and u are, respectively, the density, 
pressure and specific thermal energy of the gas; v g is the fluid 
velocity; $ is the gravitational potential; 7 is the adiabatic 
index (assumed to be 5/3 in the following); F and A are the 
specific heating and cooling rates. For spherical symmetry, 
the continuity equation and momentum equation can be re- 
written as 
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where M(< r g ) is the total mass (including dark matter) 
enclosed by a sphere with radius r g and G is gravitational 
constant. Note that the specific thermal energy of the gas 
can change not only by cooling and heating but also by 
adiabatic expansion and compression. Instead of the energy 
equation, one may choose to integrate the entropy equation. 
We define the specific entropy as 
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where T is the temperature, and n is the particle number 
density of the gas. We can then write the equation of state 
as 



where the adiabat A is related to S as 
k 



(8) 



(9) 



where k is the Boltzmann constant, /i is the mean molecular 
weight and rrin is the mass of a hydrogen atom. The adiabat 
can be changed when the gas is heated or cools, or as it is 
shocked. We can therefore write 



dA = 



(7 - l)dQ 



7-1 

Pg 



AQ = dW B hock + dQhc , 



(10) 
(11) 



where dWBhock is the work done by shock, and dQhc is the 
heat due to other cooling and heating processes. 

In most of the simulations presented in this paper, 
atomic radiative cooling is incorporated using the cooling 
function given in Katz, Weinberg & Hernquist (1996). We 
assume the gas has a primordial abundance, i.e. 76% of the 
mass is in hydrogen, 24% is in helium, and there are no heav- 
ier elements. The atomic cooling processes considered here 
work only for temperature above ~ 10 4 K. We calculate the 
abundance of different ionic species as a function of density 
and temperature assuming ionization equilibrium. Since the 
shock front is thin compared to all the length scales of in- 
terest, the mass of gas involved in the shock front is only 
a small fraction. Following Thoul and Weinberg (1995), we 
ignore cooling in the shock front. 

2.2 Numerical scheme 

We use the Lagrangian finite-difference scheme of Bow- 
ers & Wilson (1991) to follow the evolution of gas shells. 
In this scheme, phase space quantities, such as positions 
and velocities are zone-edge-centered, while thermodynamic 
quantities, such as pressure, density and entropy /energy are 
zone-centered. In what follows, subscripts denote the in- 
dex/position of gas shells. The acceleration of the ith gas 
shell is 



di = — 47rr. 



2Pi + - Pi- 
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+ 9i- 



(12) 



where p i+ and pi_ are the pressures of the gas at the outer 
and inner sides of the shell, rrii is mass of the gas shell, and 
gi is gravitational acceleration. 

The time integration is carried out by the DKD (drift- 
kick-drift) leap-frog algorithm: 

(13) 
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in which superscripts denote the current time step. Once the 
integration is carried out to the nth step, we can calculate 
the gas density between every two neighboring shells: 

3mi 



4n [ (r „ )3 „ (r n_ i)3 ] 



(16) 



During this time step, the change of adiabat is 
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where the work done by the shock is given by 
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Here p s h oc k is equivalent to the shock pressure using an arti- 
ficial viscosity technique (Richtmyer & Morton 1967; Thoul 
& Weinberg 1995): 
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< 0, and Pshock.i- = otherwise. Fol- 



lowing Thoul & Weinberg (1995), we choose a — 4. Our 
tests show that the results are not sensitive to the value of 
this parameter. The other entropy source term comes from 
radiative cooling, 

An 

AQL,,- = -- *-At n+1/2 , (20) 



where p™_ = p^p"^ 1 is the mean density of the gas shells 
during the nth time step. The advanced adiabat at the next 
time step is 

^n + l = A n_ + AA n_ _ ^ 

Other thermal quantities at the new time step are derived 
as follows: 
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We ignore radiative cooling on the shell front, at which 
v" +1 — v"^i /2 < and Vi » C„,i-, where C s ,»_ = 

\J JPi-Jpi- is the adiabatic sound speed for the shell. 

When a gas shell cools to a temperature below 10 4 K, 
the gas in the shell is considered to be cold. Since we do not 
include any cooling processes below this temperature, the 
cold gas is assumed to retain a temperature of 10 4 K until it 
flows into the center of the halo, where we drop the cold gas 
out of the gas component. Once a gas shell is dropped out, 
we do not trace its dynamical evolution anymore. Instead, 
we put the gas in a central exponential disk with a scale- 
length r d = 2^r v (Mo, Mao & White 1998), where r v is 
the virial radius of the halo at z = 0. The cold gas disk is 
assumed to be a rigid object, and its gravity is included in 
the subsequent evolution of the other mass shells. 

Each time step of the integration is controlled to be 
smaller than any of the important timescales: the dynamical, 
the Courant, and the cooling timescales, 
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where c d , cc and c c are safety parameters. Moreover, in 
order to avoid shell crossing, we have an additional time 
step control, 



2/ 



A£ ve i = min 



{ Vi-Vi-i J 



(27) 



where c v is again a safety parameter. In practice, the time 
step is chosen to be the minimum of all the above time steps, 



dt — min{d£dyn, dtcour, dt coo i, dt ve i}. 



(28) 



The safety parameters are taken to be Cd = 0.01, cc = 
0.2, c c = 0.1, and c v = 0.05, as suggested by Thoul & Wein- 
berg (1995). 

The numerical code was tested against idealized prob- 
lems with analytical solutions, such as the Sedov solu- 
tion and the self-similarity models of Bertschinger (1983). 
The code also reproduces the results of Thoul & Weinberg 
(1995). 



2.3 Initial conditions 

The simulation starts from an initial condition, which spec- 
ifies the position, velocity and other properties of each shell 
at a high redshift, z = Zi. We set up our initial condition 
based on the mass accretion histories (MAHs) of dark mat- 
ter halos. TV-body simulations show that the mass accretion 
histories of CDM halos have some universal behaviors (e.g. 
Wechsler et al. 2002; Zhao et al. 2003a; b), and can be de- 
scribed by a simple functional form, 



M(a) = M exp 



-2a c («> 



(29) 



where a is the expansion scale factor, a c is the scale factor 
corresponding to the formation time of the halo, and Mo 
is mass of the halo at the observation time cio (Wechsler et 
al. 2002). Note that a c is the single free parameter. 

Once the mass accretion history is specified, we can use 
it to determine the density profile of the initial perturbation. 
In an expanding universe, a mass shell turns around and 
then collapses at a time that depends on the initial over- 
density within the mass shell. Based on a spherical collapse 
model, a mass shell collapses at a time when the linear over- 
density reaches ~ 1.686. Thus, the initial over-density within 
a mass shell that collapses at a redshift z can be written as 



Si(M) = 1.686 



D(zj) 
D(z)> 



(30) 



where D(z) is the linear growth factor, and M is the mass 
within the mass shell. The mass M is related to the initial 
radius n of the mass shell by 

-, 1/3 

(31) 



Ti = 



M 



\np(zi)[l + 5i(M)] 



where p; = pcrit.o ■ !)m • (1 + z;) 3 i s the mean matter density 
at redshift Zi, p C rit,o is the critical density of the universe at 
the present time, and JIm is the matter density parameter. 
For a given cosmology and halo mass accretion history, the 
above relations allow us to specify Si as a function n. 

The initial velocity of each mass shell has two compo- 
nents, the Hubble expansion Vi = ViH(zi), and the peculiar 
velocity owing to the density perturbation: 



3HQ f 



S, 



where g is the peculiar gravitational acceleration, 



g(x) = Gp b a 



J dV<J(x', 



t)(x-x')/|x'-x| : 
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and / = d In D(z)/d In (1 + z), with all the quantities eval- 
uated at the initial redshift Zi. We choose z, = 200, and 
the initial temperature of gas shells is set to be the CMB 
temperature at this epoch. 

Our simulations cover a large range of halo mass, from 
10 10 to 10 15 hT 1 Mq. In hierarchical structure formation, the 
value of a c is correlated with halo mass. In our simulations, 
we use the typical value of a c for halos of a given mass. The 
parameters of the halo models are described in Table 1 . As 
in Lu et al. (2006) , particles accreted before z c are assumed 
to have isotropic velocities, so that the final dark halos have 
the Navarro, Frenk & White profiles (Navarro et al. 1996, 
1997). 

In all the simulations, we assume Qm = 0.3 and = 
0.7 at the present time. We fix the baryon fraction to be 
fb = 0.17 for our fiducial model, but vary this parameter in 
some of our simulations. We use 5 x 10 4 shells of the same 
mass for the dark matter, and 500 equal-mass shells for the 
gas component. To avoid a numerical instability caused by 
dark matter shell-crossing, the mass of each dark matter 
shell is chosen to be smaller than that of the gas shell. Our 
tests with higher resolutions show that our resolution en- 
sures numerical convergence. 



3 ACCRETION BY DARK MATTER HALOS 
IN A PREHEATED MEDIUM 

In this section, we present simulation results that neglect 
radiative cooling. First, we consider models in which the 
intergalactic medium is cold. During the formation of a dark 
matter halo, gas is accreted and shocked. In the absence of 
radiative cooling, the temperature of the post-shock gas is 
roughly given by the in-fall kinetic energy of the gas: 



T — 



16k 



(34) 



where the shock velocity v s is approximately the virial ve- 
locity of the halo at the time in consideration. Since the 
gravitational potential of the halo is established at an early 
stage of halo formation, the virial velocity does not change 
much at later times. This is the main reason for the rather 
flat temperature profiles shown in the left panel of Figure 1. 
Note that, as the halo mass increases with time, the radius 
within which the gas is shock heated also increases. 

The right panel of Figure 1 shows the evolution of the 
specific entropy profile. The specific entropy of the accreted 
gas increases with radius, because the density of shocked 
gas is lower at larger radius. The entropy profile is roughly 
a power-law and shows a self-similar behavior. Indeed, if 
scales the radius r in units of the present virial radius, and 
scales the specific entropy in units of the "virial entropy" , 
defined as 



S v 



2/3 



(35) 
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Table 1. Halo models 
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Figure 1. The evolution of temperature (left panel) and normalized entropy profile (right panel) of a M = 10 12 h^ 1 Mq halo in the 
absence of preheating. The dotted, dashed, and solid curves correspond to z = 3, 1.5 and 0.85, respectively. The radius is scaled to the 
virial radius of the halo at z = 0. The sharp breaks correspond to the shock fronts. 



where n v is the mean particle density within the virialized 
halo, the entropy profiles for halos of different masses all 
obey: 

S(r) 



(£)"■ 



S ~\r J ' (36) 

Ov \ / v / 

with a ~ 1.3, as shown in Figure 2. This result is consistent 
with the analytical work of Tozzi & Norman (2001), but the 
slope is steeper than the observed entropy profiles of nearby 
cooling flow clusters (PifTaretti etal 2005), which have a ~ 
0.95. As we shall see, radiative cooling changes the entropy 
profile. 

Next, we consider cases where the gas is preheated to 
a finite entropy before it is accreted into dark matter ha- 
los. In the simulations, the preheating entropy is gradually 
added into gas shells at early epochs, from z = 100 to 10, 
to make the numerical calculations stable. Here again we 
neglect radiative cooling. In Figure 3 we show the entropy 
profiles for cases with different initial entropies. One can 
clearly see the imprint of the initial entropy. The entropy 
profile is fiat at small radii, with an amplitude comparable 
to the initial entropy, since the entropy generated by ac- 
cretion shocks is smaller than the initial entropy at small 
radii. Since the entropy generated by shocks increases with 



radius, the entropy profile will eventually be dominated by 
shock generation and follow that of no preheating case. This 
transition occurs roughly at the radius where the shock en- 
tropy equals to the preheating entropy. The final entropy 
profile can, therefore, roughly be described by 



S(r) = 



S°(r) 



if 
if 



r < r t ; 
r > r t , 



(37) 



where S p h is the preheating entropy, S°(r) is the entropy 
expected from cold accretion, and rt is the radius at which 
S {> = Sph- Since S (r)/S v is roughly a scale-free function 
of r/r v , the characteristic radius rt, expressed in units of 
rv, is determined by S p h/S v . Hence, the effect of preheating 
in adiabatic accretion is almost entirely determined by the 
ratio Sph/SV- 

As expected, preheating reduces the amount of gas 
that accretes into a halo. In Figure 4 we show the ratio of 
fg(S P h)/fg(0) as a function of Sph/Sv, where f g (S p h) is the 
ratio between the mass of accreted gas and the halo mass, 
and / 9 (0) is the corresponding ratio for the non-preheating 
model. As one can see, the / 9 (S p h)// 9 (0) - S p h/S v relation 
is independent of both the halo mass and of the initial bary- 
onic fraction. The relation can be well described by 
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Figure 2. The normalized specific entropy profiles for halos of 
different masses at z = 0. The solid, dashed and dotted curves are 
for halos with masses of 10 10 , 10 12 and 10 14 /i -1 Mq, respectively. 
We scale the radius to the virial radius, r v , and scale the specific 
entropy to the virial entropy, S v , as defined in the text. 




Figure 3. The entropy profiles in two preheating models: the 
dotted line corresponds to S p h = 0.3S V and the dashed line cor- 
responds to S ph = O.liSV, compared to the model without pre- 
heating (solid line). We ignore radiative cooling. 



f 9 lfM = 



i + 



1/2 



(38) 



This relation shows that the reduction in the accreted gas 
becomes important when the preheating entropy is com- 
parable to or larger than S v . For S p h ^> S v , f g /f g (0) oc 
(iSph/Sv) -3 ^ 2 , and if S p h = constant then f g / f g (0) is pro- 
portional to the halo mass, M. In the last relation we as- 
sumed the virial density n v is independent of the halo mass. 
This relation is similar to that obtained by Mo & Mao (2002) 
based on a simple analytical model. 



4 COOLING OF PREHEATED GAS IN DARK 
MATTER HALOS 

Since preheating modifies the distribution of the accreted 
gas in dark halos, it can change the gas cooling rate (Mo 
& Mao 2002). In this section, we explore how preheating 
affects gas cooling in different halos. 
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4.1 Cooling in single-phase media 

To begin, we assume that the accreted gas in a dark halo has 
a single-phase. Thus, at any given time, each gas shell has a 
unique temperature and density. In the simulation, the ac- 
creted gas dissipates its thermal energy via radiative cooling, 
and its temperature drops. We consider the gas to be in a 
cold component when its temperature drops below 1 4 K. To 
characterize the efficiency of cooling, we study the fraction 
of the accreted gas that is in the cold component at z — 0. 
For simulations with preheating, we denote the total cold 
gas mass and the total accreted gas mass by M coo i(S p h) and 



Figure 4. The gas fraction that can be accreted into a halo in 
the presence of preheating, normalized by the fraction in the ab- 
sence of preheating, as a function of normalized preheating en- 
tropy, Sph/Sv Different symbols represent simulations with var- 
ious halo masses and baryon fractions, as indicated in the panel. 
The dashed curve corresponds to equation (38), while the solid 
line corresponds to the scaling relation, oc S~^ 2 . 
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Figure 6. The critical entropy Sc, as defined in the text, versus 
halo virial entropy S v . The dashed line shows S c oc S v . 

M m:c (S p h) , respectively. The corresponding quantities in the 
non- preheating cases are denoted by M coo i(0) and M aC c(0). 
Figure 5 shows the ratio / C ooi(S p h) = M cool (S ph ) / M acc (0) 
and fa,cc(S p h) = M acc (S'ph)/Macc(0) as a function of the 
preheating entropy, S p h, for halos of different masses. As 
one can see, the behavior of / CO oi(>Sph) is quite different for 
different mass halos. In massive halos, the accreted gas frac- 
tion decreases slowly with increasing amount of preheating 
entropy, while the cold gas fraction drops rapidly. Thus, in 
these systems, a large amount of the accreted gas stays in 
the hot phase. The situation is quite different for low mass 
halos, where almost all of the accreted gas can cool even if 
the IGM is preheated. 

To quantify the importance of preheating on gas cool- 
ing, we define a "critical" specific entropy S c by 

/cool (ScV/cooi (0) = 1/2, (39) 

where / CO oi(S) is the fraction of the gas that can cool in a 
halo if the IGM is preheated to a specific entropy S. Figure 6 
shows the S c — S v relation obtained from our simulations. 
Clearly, one needs a higher initial entropy to prevent the ac- 
creted gas from cooling in more massive systems. However, 
the critical specific entropy, S c , does not increase as fast 
as SV: the logarithmic slope of the S c — S v relation is ap- 
proximately 2/3 not 1. Thus, if there were a mechanism that 
could preheat the gas associated with halos to an initial spe- 
cific entropy that was proportional to S v , the effect on gas 
cooling would be more important for more massive systems. 
On the other hand, if the IGM were everywhere preheated 
to the same level of specific entropy, then the effect would 
be more important for lower-mass systems. 

It is also interesting to look at how preheating affects 
the gas cooling rate as a function of time. Since preheat- 
ing flattens the gas density distribution, the time scale for 
gas cooling is similar over a large range of radii. This is 
different from the case without preheating, where the cool- 
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Figure 7. The gas cooling rates in a 10 12 h~ 1 Mq halo versus 
rcdshift. Different symbols represent results for different preheat- 
ing levels, as indicated in the panel. 

ing time scale is shorter in the inner, denser regions and 
the cooling proceeds in an inside-out fashion. In Figure 7, 
we show the cooling rate as a function of redshift in a 
halo of 10 12 h~ x Mq. Here we consider two models of pre- 
heating, one with S p h = 10 keV cm 2 and the other with 
Sph = 14keVcm 2 , and compare the results of these two 
models to the case with S p h = 0. Againi, we assume the 
accreted gas has a single phase. As one can see, preheat- 
ing causes a delay in gas cooling. The cooling rate peaks at 
z ~ 2 in the absence of preheating, and the cooling rate at 
the present time is about 4M©/year. Preheating decreases 
the cooling rate at high z but enhances the cooling rate at 
lower z, because much of the gas has to wait until a later 
time before it can cool. For all the simulations we analyzed, 
the delay in gas cooling is more important if the ratio S p h/S v 
is larger. 

Based on the results presented above, one can conclude 
that the accretion of cold gas into galaxies occurs rather 
late in the presence of preheating. This may be a problem, 
because it implies that most of the stars in present spiral 
galaxies must form very late. As we will see below, this 
problem can be alleviated by assuming that the gas is in 
a multi-phase rather than in a single-phase state. 

4.2 Cooling in multi-phase media 

There are many reasons to believe that the gas in galaxy ha- 
los has a multi-phase structure. First of all, at the temper- 
atures considered here, the gas is expected to be thermally 
unstable, and so a multi-phase medium can develop through 
thermal instability (Field 1965; Fall & Rees 1985; Murray & 
Lin 1990). Secondly, the density field of the CDM universe 
is expected to be clumpy, which can promote thermal insta- 
bility. The development of a multi-phase medium will affect 
the cooling of the halo g clS, clS the cooling time in the high- 
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Figure 5. Diamonds show the amount of accreted gas in the presence of preheating, normalized by the amount of accreted gas in the 
absence of preheating, versus the preheating entropy. Crosses show the amount of cold gas, normalized by the amount of accreted gas in 
the absence of preheating, versus the preheating entropy. Different panels show results for halos of different masses. 



density/low-temperature phase may become shorter than 
the age, even if the single-phase cooling time does not. In 
this section, we use simple models to demonstrate how gas 
cools in a multi-phase medium. 

To incorporate a multi-phase medium in our simula- 
tions, we split each gas shell into N phases over a range of 
gas densities. The gas in each phase is assumed to move with 
its host shell and to be in pressure equilibrium with the other 
phases in any subsequent evolution. Each sub-shell, which 
represents gas in a phase, is characterized by its number 
density mj , temperature Tij and mass rriij , where i labels a 
shell, and j labels a sub-shell phase. The thermodynamics 
of the multi-phase gas is treated in the same way as that 
described in Thomas (1988). At each time step, the radia- 
tive cooling of every sub-shell is calculated using its current 
density and temperature, and a new density and a new tem- 
perature are determined through the hydrodynamical equa- 
tions (see Thomas 1988 and its appendix for details). Once 
the temperature of a sub-shell drops below 10 4 K, this phase 
drops out of the shell, and the gas in the sub-shell is then 
deposited into a central disk that is assumed to form at the 
center of the host halo. The remaining phases in the shell 
are assumed to expand adiabatically to occupy the extra 
volume. In practice, we adopt the multi-phase model for a 
gas shell as soon as it is accreted within the current virial ra- 
dius. We, therefore, neglect the details of how a multi-phase 
medium actually develops. 

The density distribution of the multi-phase gas is de- 
scribed by the volume fraction function, /(p, p)dp, which 
describes the volume fraction occupied by a phase with a 
density in the range p±dp/2 for a gas shell with mean den- 
sity p. By definition J f dp — 1. In our simulations, each 
gas shell initially contains 50 sub-shells equally spaced in 



log p. Since the form of f(p,p) is not known a priori, we 
consider three simple models, as examples, which are plot- 
ted in Figure 8. Two of these are log- normal functions with 
two different widths, cr logp / po = 0.15 and 0.05, where p is 
the median density of the distribution, related to the mean 

— — 2 /9 

density, p, by p = pe CT ' . In the third model, the volume- 
fraction function is assumed to have the "cooling tail" form, 

/( P ,p)ocp-< 4 - Q >exp[-(p7p) 2 -"] , (40) 

where a is the slope of the cooling function, A a T°. This 
model has been used to describe multi-phase cooling flows 
in galaxy clusters (Nulsen 1986; Thomas 1988; Waxman & 
Miralda-Escude 1995). We have adopted a = 0.4, which 
roughly describes the cooling by free-free process. In the 
simulations, we only sample phases with / > 10~ 5 . Our 
test with the "cooling tail" model shows that our simulation 
reproduces the results obtained by Thomas (1988). 

To demonstrate the effect of changing the gas from a 
single-phase medium to a multi-phase medium, we consider 
a case where the final halo mass is 10 1 h~ Mq, and the ini- 
tial preheating entropy S p h = lOkeVcm , using the three 
models for / described above. Figure 9 shows the gas cooling 
histories, i.e. the mass of cooled gas as a function of redshift, 
compared to that of a single-phase model. As expected, the 
cooling history is more extended in models where the den- 
sity distribution is broader. For a broader density distribu- 
tion, there is larger amount of high-density gas that can 
cool rapidly. If the density distribution is as broad as the 
"cooling-tail" model, cold gas can be assembled into galax- 
ies much earlier than in the single-phase model. 

Remember that in the single-phase model, cooling is 
inside-out. So if the cooling time of the gas were comparable 
to the age, only gas in the innermost region could cool. On 
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Figure 8. The volume fraction occupied by phases of different 
density. The dotted and dashed curves are log-normal distribu- 
tions with dispersions cr logp ^ pg = 0.05 and 0.15, respectively. The 
solid curve is the "cooling tail" model described in the text. The 
x-axis is scaled by the mean density, p. 



the other hand, in a multi-phase medium, cooling can occur 
almost everywhere in the high-density phases. To demon- 
strate this, we plot in Figure 10 the initial (Lagrangian) radii 
that contain 20%, 50% and 80% of the gas that drops out 
to the cold phase in a time interval At = 0.5 Gyrs as a func- 
tion of (l + z). In the single-phase case, all the gas within 
a certain Lagrangian radius can cool, and so the 20%, 50% 
and 80% radii are all the same (the dashed curve in Fig- 
ure 10). For the single-phase case shown in Figure 10, this 
"cooling radius" is much smaller than the virial radius of 
the halo (dot-dashed curve), implying that only gas in the 
innermost parts of the halo can cool. In the multi-phase 
case, however, at any time gas can cool to the cold phase 
over a large range of Lagrangian radii. This result could have 
important implications for the formation of galaxy disks. N- 
body simulations show that the specific angular momentum 
of dark matter particles increases roughly linearly with the 
Lagrangian radius (Bullock 2001). Thus, if the angular mo- 
mentum of the gas follows that of dark matter, the specific 
angular momentum of the cooled gas in the single-phase case 
is expected to be lower than the average specific angular mo- 
mentum of the dark matter, since only the gas in the inner 
regions can cool. This could lead to the formation of galaxy 
disks that are too compact (e.g. Mo, Mao & White 1998). 
On the other hand, in the multi-phase model, the cold gas 
that settles into the disk comes from various radii, and so 
its average specific angular momentum may be comparable 
to that of the dark matter, a condition that is required to 
reproduce the observed sizes of galaxy disks (e.g. Mo, Mao 
& White 1998). 

As demonstrated in Figure 11, if the initial specific en- 
tropy, /Sph, is comparable to or smaller than the critical 
value, 5 C , the total amount of gas that can cool at the 




1 10 

l+z 

Figure 9. The total amount of cooled gas in a 10 12 Mq 
halo with preheating entropy S p h = S c ~ lOkeVcm 2 . The dash- 
dotted line shows the result of a model that assumes the halo 
gas is a single-phase medium. The three other curves show the 
results of the three multi-phase models: the cooling-tail model 
(solid curve), the log- normal model with a = 0.15 (dashed curve), 
and log-normal model with a = 0.05 (dotted curve). 

present time is quite insensitive to the assumed gas den- 
sity distribution, although different models of the gas den- 
sity distribution produce different accretion histories of cold 
gas. This occurs because the cooling time of the gas accreted 
into halos, even when it is assumed to be in a single-phase, 
is comparable to or shorter than the age of the universe 
at the present time. However, if the value of S p h is much 
larger than S c , the total amount of gas that cools by the 
present time can be significantly larger in the multi-phase 
model than in the single-phase model, for halos with masses 
>1O 11 /i _1 M . For halos with smaller masses, cooling is so 
effective that almost all the accreted gas can cool even if 



5 APPLICATIONS 

So far we have discussed in general how preheating af- 
fects gas accretion and cooling in dark matter halos. In the 
present section, we consider two specific examples that may 
be relevant to the actual formation of galaxies. 

5.1 Preheating by a central AGN 

First, we consider the case where the gas associated with a 
dark halo is heated by the outflow of a central AGN. Re- 
cent observations show that almost all galaxies, at least 
those with a significant spheroidal component, i.e. ellip- 
tical galaxies and spiral galaxies with a significant bulge, 
host a central supermassive black hole (SBH) with a mass 
that is correlated with the velocity dispersion of the bulge, 
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Figure 11. The mass of cooled gas as a function of redshift. In each panel, the solid curve shows cooled gas mass in a simulation assuming 
a multi-phase medium with the "cooling tail" density distribution, and the dashed curve stands for the result of the corresponding single- 
phase simulation. The dash-dotted curve shows the baryon fraction (/;, = 0.17) times the halo mass. Results are shown for halos of 
different present masses which are noted in the right column panels for each row and different preheating entropies which are noted at 
the bottom of each panel. 



M BH = 10 s - 3 M Q (cr/200kms- 1 ) 4 02 (Tremaine etal 2002). 
This correlation is quite independent of redshift out to z ^ 3 
(Sheilds et al. 2003). If fed with gas, these black holes are 
able to power AGNs that may produce outflows and heat the 
gas in dark matter halos. This kind of feedback process has 
been suggested as a plausible mechanism to quench both 
star formation and black halo growth in massive galaxies 
(Binney & Tabor 1995; Silk & Rees 1998; Nulsen & Fabian 
2000; Springel et al. 2005). However, it is still unclear how 



the energy associated with the AGNs is coupled to the gas, 
allowing it to quench gas cooling and star formation. It has 
been suggested that the energy deposition is localized so 
that the outflows driven by AGNs are thermalized in the 
ambient medium only within the host galaxy (e.g. Springel 
et al. 2005). However, recent numerical simulations with jets 
show that the jet may tunnel through low-density regions of 
the ambient medium, and that the energy may be easily 
carried to large distances and thermalized there (Vernaleo 
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Figure 10. The Lagrangian radius of the gas shell that cools as 
a function of redshift normalized to the present time (z = 0). 
Results are shown for a 10 12 /i -1 Mq halo with preheating en- 
tropy Sph = S c = lOkcVcm 2 , and the "cooling-tail" model is 
adopted for the multi-phase medium. The solid curve shows the 
Lagrangian radius within which 50% of the cooled gas is con- 
tained at each redshift, while the two dotted curves mark the radii 
which contain 20% and 80% of the cooled gas. For comparison, 
the dashed curve shows a model in which the medium is assumed 
to be single-phase. The dot-dashed curve is the Lagrangian radius 
of the forming halo. 

& Reynolds 2005). Without going into details, we consider 
a simple model here. 

We assume that the energy injected by an AGN is pro- 
portional to the mass of the SBH, £agn oc Mbh- Obser- 
vations of spiral galaxies with SBHs show that the mass of 
the SBH is tightly correlated with the circular velocity of the 
galaxy, Mbh oc (Ferrarese 2002). Since v c is related to the 
virial mass of the spiral galaxy's host halo, M oc v% , we have 
Eagn/M oc t> 2 . Thus, if the energy of the AGN is to heat 
all the gas associated with its host halo, the energy gained 
by each particle from the AGN feedback is roughly propor- 
tional to « 2 , i.e. to the virial temperature of the host halo. If 
such feedback occurs at the same redshift for systems with 
different masses, the specific entropy owing to AGN preheat- 
ing is expected to be proportional to the virial entropy of 
the halo, i.e. S p h oc S v (Scannapieco & Oh 2004). However, 
observations of quasars suggest that there is an upper limit 
to quasar luminosities (Fan et al. 2003), corresponding to 
a SMB mass of Mbh ~ 10 9 M Q . This mass corresponds to 
v c = 500km s -1 assuming Eddington luminosities and the 
observed Mbh — v c relation (Wyithe & Loeb 2003). Such 
a limit in Mbh implies that the Mbh-v c relation must be 
broken for massive halos. There must be some process that 
quenches the growth of the SMB and, given the Mbh-t rela- 
tion, also the growth of the central galaxy in massive halos. 
In what follows, we will consider a preheating model with 
Sph oc S v , and examine if the observed mass limit can be 
accommodated. 



Based on the above discussion, we write 
Sph — ff p h,i2 „ V , (41) 

<-5v,12 

where S v ,i2 is the virial entropy of a 10 12 h^ 1 Mq halo and 
<Sph,i2 is a normalization factor, taken to be the value of 
Sph corresponding to that in a halo of present mass M = 
10 12 h~ x M©. To estimate the normalization factor, we use 
the model presented in Scannapieco & Oh (2004). Since the 
lifetime of an AGN is typically ~ 10 7 years, much shorter 
than a Hubble time, the effect of an AGN on gas at large 
distances may be modeled as a blast wave. Adopting the 
Sedov- Taylor solution, we can write the postshock entropy 
as 

S s = 1.8 x 10 4 keVcm 2 £ 60 57 5/3 (l + z)~ 6 R~m pc , (42) 

where Eeo is the kinetic energy in an AGN outflow, E, in 
units of 10 60 ergs, 5 B is the overdensity of the surrounding 
medium defined as 8 B — pb/pb, and i? s ,Mpc is the shock radius 
in Mpc and is given by the Sedov solution as 

R B = 1.7Mpc^ B *r 1/B (l + *r S/5 C> («) 

with tc yT , the evolution time of the shock in gigayears. 
Following Scannapieco & Oh (2004), we assume that E 
is a constant fraction of the bolometric luminosity, Lboi: 
E ~ ek-^boiiAGN, where we take £k ~ 0.05. The AGN 
life time, £agn, is assumed to be proportional to the dy- 
namical time of a virialized halo, and is approximated by 
Ugn ~ 5.2 x 10 7 yr (1 + z)~ 3/2 (e.g. Barkana & Loeb 2001). 
We assume the bolometric luminosity related to the lumi- 
nosity in the blue band, Lboi = 10.4Lb, as suggested by 
observations (Elvis et al. 1994). Assuming that an AGN 
shines at its Eddington luminosity, one can relate the B- 
band luminosity to the mass of the black hole (Mbh) as 
^Edd.B = 5.73 x 1O 3 L Mq 1 Mbh (e.g. Willot, McLure & 
Jarvis 2003) . We use the observed relation between Mbh and 
the circular velocity v c obtained by Ferrarese et al. (2002), 

Mbh = T4x1O«M f(^- t ) 5 , (44) 

to estimate the black hole mass, where F is a model param- 
eter taken to be 0.6, as in Wyithe & Loeb (2003). With all 
these assumptions, one can estimate Sph, 12 ~ 7keVcm 2 at 
2 = 2, assuming 5 S = 1. 

We first carry out single-phase simulations similar to 
those presented in Section 4, but with the preheating en- 
tropy set by (41) with Sph, 12 = 7keVcm 2 . For compar- 
ison, we also consider another sequence with S p h,i2 = 
3.5 keV cm 2 . Figure 12 shows the mass of cooled gas at 2 = 
as a function of halo mass. It is clear that the suppres- 
sion of gas cooling by AGN preheating is more important 
in more massive halos. This result can be understood us- 
ing Figure 6. Since the AGN preheating entropy considered 
here is proportional to the virial entropy of the halo, the 
effect of preheating is more important in more massive ha- 
los. In this model, there is an upper-mass limit for the 
cooled gas mass, which is about ~ 3 x 10 11 hT 1 Mq assuming 
Sph, 12 = 7keVcm 2 , and about ~2x 10 12 hT 1 M Q assum- 
ing Sph, 12 = 3.5keVcm 2 . This maximum mass occurs at a 
halo mass of M ~ 4 x 10 12 for S p h,i2 = 7keVcm 2 , and at 
M~5x 10 13 for S p h,i2 = 3.5keVcm 2 . Thus, for a given 
preheating amplitude S p h,i2, significant gas cooling in halos 
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above a certain mass can only occur in progenitors before 
the formation of the central AGN. As soon as the central 
AGN has formed, gas cooling is effectively suppressed by 
preheating. 

This result has significant implications. Observations of 
AGNs suggest that there is an upper limit for quasar lumi- 
nosity (Fan et al. 2003), which corresponds to a black hole 
with a mass of ~ 10 9 Mq assuming it radiates with an Ed- 
dington luminosity. This mass is approximately equal to that 
of the central black hole in a halo with a circular velocity 
v c ~ 470km s _1 (corresponding to a halo mass of 10 13 - 5 M 
at the present time), using the observed Mbh — v c relation 
(Wyithe & Loeb 2003). If the mass of a central black hole 
scales with its host halo as the Mbh — v c relation for all dark 
matter halos, such an upper limit for black hole mass would 
not be expected. Thus, there must be a process that can pre- 
vent the growth of central black holes in massive halos. The 
AGN preheating model considered here may alleviate this 
difficulty. Before the formation of the supermassive black 
hole, gas can cool effectively and the mass of the central 
black hole grows rapidly. However, as soon as the mass of the 
central black hole reaches a value of ~ I0 9 M Q , gas cooling is 
effectively suppressed and the black hole stops growing. The 
same process also prevents the growth of the central galaxy, 
so that the 1-d velocity dispersion of the central galaxy can- 
not be much larger than 330km s -1 . This may explain why 
not many galaxies can have velocity dispersions much larger 
than this (e.g. Sheth et al. 2003). 

There are, however, a number of uncertainties in this 
conclusion. First of all, if the growth of black hole mass is 
indeed saturated at a fixed value, the preheating produced 
by the AGN would become unimportant as the dark matter 
halo grows and the corresponding S v becomes larger than 
the entropy generated by the AGN. Cooling will then resume 
in halos with v c S> 330km s -1 . To prevent further cooling in 
massive halos after the formation of SMBs, one may still 
need energy sources to balance the radiative cooling. One 
possibility here is the feedback from the central galaxy in 
"radio" mode. Although the energy feedback in such a mode 
is much weaker than that in the "quasar" mode, it may be 
sufficient to balance radiative cooling in massive halos (e.g. 
Croton et al. 2006). Another uncertainty is the possible de- 
velopment of a multi-phase medium in massive halos, which 
may promote additional radiative cooling. To demonstrate 
this, we run multi-phase simulations, assuming the "cool- 
ing tail" model, for halos more massive than 10 12 hT 1 Mq. 
The dashed curves in Figure 12 show the mass of cooled gas 
at 2 = as a function of halo mass for this model. Now, 
there is no sharp truncation in the cold gas mass, although 
the total cold gas mass is still reduced relative to the non- 
preheating model. Unfortunately, it is unclear whether a sig- 
nificant multi-phase medium can develop in massive halos. 
Indeed, in order for the 'radio-mode' feedback to work in 
massive halos, the formation of multi-phase media has to be 
suppressed. 

5.2 Low-mass halos in a warm intergalactic 
medium 

As one can see from Figure 12, AGN preheating has little ef- 
fect on small halos, because the black holes in these systems 
are small and the cooling is efficient. Thus, AGN preheating 




M (h _1 M.) 

Figure 12. Cooled gas mass versus halo mass in simulations 
assuming AGN preheating. The squares show results assuming 
<S'ph,l2 = 7keVcm 2 , and the triangles show results assuming 
S'ph,l2 = 3.5keVcm 2 . Symbols connected by solid lines are re- 
sults assuming a single-phase medium, and symbols connected 
by dashed curves arc results assuming the "cooling tail" model 
of the multi-phase medium. For comparison, diamonds connected 
by dotted lines are the results assuming no preheating. 



cannot suppress star formation in low-mass halos to explain 
the flat faint-end slope of the galaxy luminosity function. 
Some other processes must be at work in small halos. In 
this subsection, we consider a model in which all small ha- 
los (M < 10 12 M Q ) are assumed to be embedded in a warm 
medium with an initial specific entropy S p h ~ lOkeVcm 2 . 
This model is motivated by the proposal of Mo et al. (2005) 
that low-mass halos at low redshift may be surrounded by a 
warm intergalactic medium produced by gravitational pan- 
caking. Another possibility for such preheating is supernova 
explosions associated with the formation of halo stars, bulge 
stars, and stars in the thick disk. 

We consider a case where S p h = 8 keV cm 2 . Since such a 
low value of S p h has little effect on halos with M > 10 12 Me, 
we concentrate on smaller halos. Figure 13 shows the mass of 
cooled gas, M coo i, as a function of halo mass. For low mass 
halos, M coo i scales with halo mass roug hly as M 2 . This can 
be understood as follows. As we have seen in Section 3, if 
Sph > S v , as is the case for low-mass halos, the fraction 
of accreted gas scales as / acc = (Sph/Sv) -3 ^ 2 - Since at a 
fixed redshift S v oc T v and T v oc M 2/3 , the mass of the 
accreted gas M acc = f^ccftM oc M 2 for a constant preheat- 
ing entropy. In small halos almost all the accreted gas can 
cool, so we have M coo i oc M 2 . This result supports the as- 
sumption made by Mo et al. (2005) about the M coo i - M 
relation for low-mass halos accreting gas in a medium pre- 
heated by the formation of pancakes. As shown in Mo et 
al. (2005), such a M coo i - M relation is required to explain 
bothe the faint-end slope of the galaxy luminosity function 
and the Hi-mass function of galaxies. It should be noted that 
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Figure 13. The mass of cooled gas versus halo mass for a model 
in which the intergalactic gas is preheated to a specific entropy 
S p h = 8.0keVcm 2 (stars). In this case, the cooled gas mass scales 
with halo mass roughly as M 2 . For comparison, the diamonds 
show the results for a model without preheating. 

the multi-phase model does not change the cold gas mass in 
these halos. As we have demonstrated in Section 4, the cool- 
ing is always efficient in small halos. Hence, once the gas is 
accreted into low-mass halos, it can cool no matter it was in 
high-density or low-density phases. Note, however, that the 
assembly history of the cold gas into the central disk may 
depend significantly on whether the medium is single-phase 
or multi-phase, as discussed in Section 4.2. 



6 SUMMARY 

In this paper, we examine the effects of preheating on gas 
cooling and accretion into dark matter halos, using a 1-D 
hydrodynamical code and realistic halo formation histories. 
We explore the behavior of the preheated gas by varying the 
virial mass of the halo and the amount of preheating. We 
find that preheating can both prevent gas accretion into dark 
matter halos and reduce the radiative cooling of that gas. 
If the preheating entropy is larger than the virial entropy 
of the halo, the fraction of the gas that can be accreted 
into a halo scales roughly as f g oc (S , p h/5 V ) _3/ ' 2 . For low- 
mass halos, almost all the accreted gas can cool, and so 
preheating reduces the total amount of cold gas that can be 
assembled into the central galaxy only by reducing the total 
amount of gas that can be accreted by the halo. For massive 
halos, on the other hand, preheating not only reduces the 
amount of gas that can be accreted, but also suppresses 
radiative cooling. We also examine the effects of a multi- 
phase medium on the cooling of halo gas and find that the 
cooling history of multi-phase gas is very different from that 
of a single-phase medium model. In a multi-phase medium, 
gas cooling can occur in the dense phases soon after the 



gas is accreted into the halo, in contrast to a single-phase 
medium where gas cooling is delayed. Also, in a single-phase 
medium, gas cooling in general proceeds in an inside-out 
fashion, while in a multi-phase medium gas can cool, at a 
given time, over a wide range of radii. 

We consider two examples where preheating may have 
played important roles in galaxy formation and evolution. 
First, we examine a case where the gas associated with a 
dark matter halo is preheated by an AGN outflow powered 
by a central supermassive black hole. Using assumptions 
made in recent quasar feedback models, we demonstrate that 
such preheating can effectively suppress gas cooling in halos 
with masses above ~ lO 13 /i _1 M0, leading to an upper mass 
limit of about 10 12 Mq for the central galaxy that can 
form in such a massive halo. This may be the reason why 
the galaxy luminosity (mass) function has a sharp break at 
the bright (massive) end. This suppression of gas cooling 
in massive halos may also be responsible for the observed 
bimodality in the colour-magnitude diagram. Observations 
show that the blue sequence, which consists of galaxies with 
significant star formation, is truncated at the bright end, 
with very few blue galaxies having stellar masses exceed- 
ing W 11 Mq. The more massive galaxies are all red, with 
little current star formation. This suggests that there must 
be some process to quench star formation in these massive 
galaxies. Recent semi-analytical models of galaxy formation 
including AGN feedback show that the bimodality in the 
colour distribution and the truncation of the blue sequence 
can be reproduced if the feedback from the central galaxy is 
in the "radio mode" (e.g. Croton et al 2005; Cattaneo et al 
2006). However, in order for the "radio-mode" feedback to 
be effective, the halo gas must be in a diffuse hot form. In 
the scenario we consider here, the halo gas is preheated by 
the feedback from the central AGN in the "quasar mode", 
and any subsequent cooling is suppressed in the preheated 
medium. This may help create the hot diffuse halos within 
which the "radio mode" feedback can act effectively. 

The second example we consider is the case where the 
intergalactic medium is assumed to be preheated to a con- 
stant entropy level. This model is motivated by the scenario 
outlined in Mo et al. (2005) where such preheating may help 
explain the observed stellar- and Hi-mass functions of galax- 
ies at the low-mass end in current CDM models. Our sim- 
ulations confirm that such preheating can indeed produce 
the cold-gas-mass/halo-mass relation, M coo \ oc M 2 , as is re- 
quired to match the observed stellar- and Hi-mass functions 
in the current CDM theory. 

The results obtained in this paper demonstrate that 
preheating may have played an important role in galaxy 
formation. To explore the observational consequences of the 
preheating scenario, one needs to combine the results ob- 
tained here with the dark matter halo population predicted 
by current CDM models of structure formation, to quantify 
how preheating affects the properties of the galaxy popula- 
tion and of the intergalactic medium. The best approach is 
perhaps through semi-analytic modeling, and we plan to do 
this in a future paper. 
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